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TABLE 2 

Parameters for Bisector Regressions 
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TABLE 3 
Surface Density Regression Slopes 
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TABLE 1 
Sources of Maps through 1994 
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Single beam maps: 




Single beam maps: 




Bajaja, Huchtmeier & Klein (1994) 


Effelsberg 


Briggs (1982) 


Arecibo 


Briggs et al. (1980) 


Arecibo 


Carignan (1985) 


Parkes 


Corbelli, Schneider & Salpeter (1989) 


Arecibo 


Giovanardi & Salpeter (1985) 


Arecibo 


Giovanelli, Williams & Haynes (1991) 


Arecibo 


Haynes (1981) 


Arecibo 


Haynes, Giovanelli & Roberts (1979) 


Arecibo 


Helou, Salpeter & Terzian (1982) 


Arecibo 


Helou, Hoffman & Salpeter (1984) 


Arecibo 


Hewitt, Haynes & Giovanelli (1983) 


Arecibo 


Hoffman et al. (1989) 


Arecibo 


Huchtmeier (1979) 


Effelsberg 


Huchtmeier & Seiradakis (1985) 


Effelsberg 


Huchtmeier, Seiradakis & Materne (1980) 


Effelsberg 


Huchtmeier, Seiradakis & Materne (1981) 


Effelsberg 


Huchtmeier & Witzel (1979) 


Effelsberg 


Krumm & Burstein (1984) 
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van Zee, Haynes & Giovanelli (1995) 


Arecibo 
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Begeman (1989) 


WSRT 


Bosma, van der Hulst & Athanassoula (1988) 


WSRT 


Bottema (1989) 


WSRT 


Bottema, Shostak & van der Kruit (1986) 


WSRT 


Boulanger & Viallefond (1992) 


WSRT & DAO 


Braine, Combes & van Driel (1993) 


WSRT 


Brinks & Klein (1988) 


VLA B 


Broeils (1992) 


WSRT 


Broeils & van Woerden (1994) 


WSRT 


Carignan & Beaulieu (1989) 


VLA D 


Carignan, Beaulieu & Freeman (1990) 


VLA C 


Carignan et al. (1990) 


WSRT 


Carignan & Freeman (1988) 


VLA D 


Carignan & Puche (I990a,b) 


VLA CD 


Carignan, Sancisi & van Albada (1988) 


WSRT 


Carilli & van Gorkom (1992) 


VLA C&D, WSRT 


Casertano & van Gorkom (1991) 


VLA D 


Cayatte et al. (1990, 1994) 


VLA D 


Chengalur, Salpeter & Terzian (1994) 


VLA C 


Comte, Lequeux & Viallefond (1985) 


VLA D 


Cote, Carignan & Sancisi (1991) 


WSRT 


Deul & van der Hulst (1986) 


WSRT 


Dickey, Hanson & Helou (1990) 


VLA D 


England (1989) 


VLA BC&CD 


England, Gottesman & Hunter (1990) 


VLA C 


Guhathakurta et al. (1988) 


VLA D 


Hummel, Dettmar & Wielebinski (1986) 


VLA CD 


Hunter, van Woerden & Gallagher (1994) 


VLA C&D 


Irwin & Seaquist (1991) 


VLA BCD 


Irwin et al. (1987) 


VLA C&D 


Israel & van Driel (1990) 


WSRT 


Jobin & Carignan (1990) 


VLA C&D 


Kamphuis & Briggs (1992) 


VLA D 


Koribalski et al. (1993) 


VLA CD 


Lake, Schommer & van Gorkom (1990) 


VLA C 


Lake & Skillman (1989) 


VLA CD 


Lequeux & Viallefond (1980) 


WSRT 


Li & Seaquist (1994) 


VLA C 


Liszt (1992) 


VLA D 


Lo, Sargent & Young (1993) 


VLA C&D 


Martimbeau, Carignan & Roy (1994) 


WSRT 


Ondrcchen & van der Hulst (I989a,b) 


VLA C 


Pedlar et al. (1992) 


VLA C 


Peters et al. (1994) 


TEST 


Puche & Carignan (1991) 


VLA CD & D 


Puche, Carignan & Bosma (1990) 


VLA CD 


Puche, Carignan & van Gorkom ( 1991) 
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Puche, Carignan & wainscoat (1991) 
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Fucne et al. (1992) 


VLA BtU 


Hand (1994) 


WbKl 


Rownd, Dickey & Helou (1994) 


VLA D 


Rupen (1991) 


VLA BCD 


Saikia et al. (1990) 


VLA BC 


Sandage & Fomalont (1993) 


VLA CD 


Sargent, Sancisi & Lo (1983) 


WSRT 


Schneider & Corbelli (1993) 


VLA BCD 


Shostak & Skillman (1989) 


WSRT 


Skillman & Bothun (1986) 


VLA D 


Skillman et al. (1987) 


VLA C 


Skillman et al. (1988) 


VLA D 


Tacconi & Young (1986) 


VLA D 


Taylor, Brinks & Skillman (1993) 


VLA D 


Taylor et al. (1994) 


VLA C 


Tully et al. (1978) 


WSRT 


van Albada et al. (1985) 


WSRT 


van Driel & Buta (1991) 


WSRT 


van Driel & van Woerden (1994) 


WSRT 


van Moorsel (1988) 


VLA D 
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ABSTRACT 

Several measures of galaxy size and mass obtained from the neutral 
hydrogen mapping of 70 dwarf irregular galaxies (Sm, Im and BCD, hereafter 
generically called "dwarfs" irrespective of size or luminosity) presented in 
the preceding paper are compared statistically to those for the set of all 
available H I-mapped dwarfs and H I-mapped spirals distributed within the 
same spatial volume to investigate variations in Tully-Fisher relations and in 
surface densities as functions of galaxy size and luminosity or mass. Some 
ambiguities due to the "non-commutativity" of the correlations among the 
variables are addressed. From linear regressions of logarithms we find that 
t oc r 2 ' 68 oc v 3 ' 73 oc m]f 5 oc m^y® where I is blue luminosity, r is the geometric 
mean of the radius to the outermost detectable H I and the optical radius, v 
is the velocity profile half-width incorporating rotation and random motions, 
ran is the mass of H I, and md yn — v 2 r/G. All are normalized by the values 
appropriate to a galaxy with blue luminosity of 10 9 Lq, typical of the region of 
overlap between dwarf and spiral galaxies. The surface density £# of neutral 
hydrogen (averaged within the "isophotal" radius r) is almost constant along the 
sequence of size/mass/luminosity while surface density Ex of blue luminosity 
increases with galaxy size. 

For quantities not involving H I we find no evidence for a "break" between 
dwarfs and spirals, but we do find some curvature in logv vs. logr and in the 
Tully-Fisher relation, logf vs. log£. Two consequences are: (i) I oc v 3 ' 7 is more 
appropriate for the whole sequence than is £ oc v 2 ' 5 as found for large spirals 
alone; (ii) the surface density of total dynamic mass (oc v 2 /r) is almost 
constant along the lower portion of the luminosity sequence, but increases 
appreciably with i along the upper portion. 
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There is an indication for a difference in the correlations involving H I mass 
or radius between dwarfs alone and spirals alone, in the sense that irregulars 
have somewhat more H I mass or slightly larger H I radii than spirals at a given 
blue luminosity, optical radius, or velocity profile width. It is not clear if this is 
a true morphological effect or merely due to mn varying less strongly than £. 

Subject headings: Galaxies: Irregular — Galaxies: Kinematics and Dynamics - 
Galaxies: Structure — Radio Lines: Galaxies 



1. Introduction 

In Paper I (Hoffman et al. 1996) we have presented H I mapping of a sample 
of 70 dwarf irregular galaxies distributed throughout the Virgo Cluster and the Local 
Supercluster. Here and throughout this paper, we use "dwarf" to mean morphological type 
Sdm, Sm, Im or BCD irrespective of the intrinsic luminosity or size of the galaxy. Other 
research groups have mapped dwarf galaxies in the meantime, and the available literature 
on mapped spiral galaxies has also been growing steadily. In this paper we combine 
our sample of mapped dwarfs with those mapped in H I by other authors using either 
single beam instruments or synthesis arrays, and with a sample of mapped spiral galaxies 
distributed within the same volume of space. The combined sample is of interest for several 
questions: Roberts & Haynes (1994) present a thorough study of how physical parameters 
vary along the morphological sequence. Are dwarf irregulars, statistically speaking, just the 
low mass continuation of the spiral sequence, or a population with a distinct evolutionary 
history (Hodge 1989 and references therein; Gallagher & Hunter 1984; Fanelli, O'Connell 
& Thuan 1988; Tyson & Scalo 1988; Dufour & Hester 1990; Schombert et al. 1990; 
Westerlund 1990; Drinkwater & Hardy 1991; Kruger & Fritze-von Alvensleben 1994; van 
den Bergh 1994)? Does the evidently episodic star formation history (Izotov, Thuan & 
Lipovetsky 1994; McGaugh 1994; Meurer, Mackie & Carignan 1994) of the dwarfs cause 
the Tully-Fisher correlations of optical properties (luminosity or diameter) vs. rotation 
velocity to deviate systematically from that for spirals (Hoffman, Helou & Salpeter 1988; 
Schneider et al. 1990; Gavazzi 1993; Sprayberry et al. 1995; Zwaan et al. 1995)? Do 
the unusual systems DDO 154, DDO 137 (Paper I) and H I 1225+01 (Giovanelli, Williams 
& Haynes 1991) continue to appear as unique systems when compared to a larger sample? 

Recently a sequence of Low Surface Brightness (LSB) galaxies ranging in size from 
Malin I , comparable in size to the largest spirals but with much reduced star formation 
(past and current) across the disk, down to LSB irregulars of very small size, has been 
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identified (Bothun et al. 1987, 1990; Impey & Bothun ; Sprayberry et al. 1993; 
McGaugh, Schombert & Bothun 1995; McGaugh 1995). Another recent study (van Zee, 
Haynes & Giovanelli 1995) observed galaxies with particularly low blue luminosity to H I 
mass ratios, which favors LSB galaxies. Our Paper I selected dwarf galaxies by morphology, 
but does not necessarily favor LSB. The combined sample of the present paper therefore 
covers a very large range in each of five extensive variables (blue luminosity Lb, radius R, 
velocity profile half-width V c , hydrogen mass M H , and indicative dynamic mass M dyn ) and 
of surface brightness S L = L B /^R 2 - The selection criteria for our sample are given in 
detail in Sect. 2 (and selection biases are discussed in Sect. 3.1), but some salient features 
(besides the broad range covered) are: H I mapping was required for all dwarfs and spirals. 
In addition to dwarfs (Sdm, Sm, Im and BCD), only late-type spirals (Sb through Sd) 
distributed within the same spatial volume were considered. The Virgo cluster was included 
in the volume, but the omission of SO, Sa and Sab types and the insistence on a measured 
H I radius (which selects against gas-stripped galaxies) means that the Virgo cluster core 
does not dominate our sample. 

We selected a large and representative sample, but did not insist on high precision 
measurements. Especially for the dwarfs, intrinsic variations as well as measuring errors 
are large, and some pairs of quantities are poorly correlated. We therefore pay special 
attention to the ambiguities in the correlation statistics, following the discussion of Isobe et 
al. (1990 — hereafter IFAB). Various measures of galaxy size and mass are available, and 
we select a particular set of five for analysis and discuss correlations among them in Sect. 
4. We address questions of ambiguity in Sect. 4.1 and questions of whether or not there are 
breaks between the dwarf and spiral sequences or curvature on log-log plots in Sect. 4.2. 
We shall see that some of the ambiguity arises from "non-commutativity" in correlation 
statistics and from similar effects for products and ratios: If some definition for regression 
lines gives b oc ar ', c oc a 7 , (be) oc a p , then p can sometimes differ appreciably from (/3 x 7). 
Our data do not allow us to obtain local surface density profiles, but we have El, E# 
and T,dyn, the surface densities of blue luminosity, H I mass and indicative gravitational 
mass averaged over the entire disks of the galaxies. The correlations among these three 
"intensive" properties, and their variations with the extensive variables, are explored in 
Sect. 5. The variation of the ratio E^/E^n = Ls/Mdyn is particularly interesting and 
particularly controversial. We compare our analysis with that from four previous surveys 
(Gavazzi 1993; Roberts & Haynes 1994; Sprayberry et al. 1995; Zwaan et al. 1995) and 
discuss implications in Sect. 6. 
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2. Sample Selection 

The definition of the sample we mapped at Arecibo is given in Paper I. In brief, the 
galaxies are chosen from three lots: 45 from the Virgo Cluster Catalog (Binggeli, Sandage & 
Tammann 1985) and Leo (Ferguson & Sandage 1990); 14 from the field survey of Binggeli, 
Tarenghi & Sandage (1990); and an additional 12 chosen to complete the sample of nearby- 
dwarfs within 6 Mpc of the Sun mapped at Arecibo. In combination, these subsamples 
produce a rather heterogeneous lot. 

Additional dwarfs mapped in H I were culled from the available literature (through 
the end of 1994). All dwarfs mapped with sufficient resolution to determine H I extents 
have been included, whether from single beam or synthesis array mapping. We adopted for 
a sample of late-type spirals (types Sb through Sd) all those mapped similarly within 20 
Mpc of the Local Group (H Q = 75 km/s/Mpc), excluding those thought to be members of 
tidally interacting pairs or groups. Neither sample is complete in any sense (except that we 
have not rejected any non- interacting mapped galaxy a priori). There is some bias toward 
systems in the Virgo cluster and its environs since so much Arecibo mapping effort has been 
expended in that area of the sky (Helou et al. 1981; Helou, Hoffman & Salpeter 1984; 
Hoffman et al. 1989), but the majority of both samples (dwarfs and spirals) is distributed 
throughout the Local Supercluster. In all, including our 70 galaxies, we find a total of 
112 mapped dwarfs out to a distance of 20 Mpc from the Local Group with another 16 at 
greater distances; the sample of mapped spirals numbers 119. We drew the data from the 
original references (see Table 1) in each case. 

In all cases, we have sought to obtain measures of galaxy properties equivalent to those 
for our Arecibo mappings as defined in Paper I: blue luminosity L B , optical radius at the 
25 mag arcsec -2 isophote i?25, dynamical speed (rotation plus dispersion) V c , galaxy-wide 
H I mass Mh, H I radius at 1/e of the central H I intensity Ru,e and H I radius at the 
outermost detected point (at a sensitivity of typically a few xlO 19 atoms cm -2 ) Rn,max- 
Optical properties were taken from RC3 (de Vaucouleurs et al. 1991) or equivalent 
photometry presented in the literature for each galaxy. For the dynamically-relevant 
velocity we adopt V c 2 = V rot 2 + 3a z 2 where V rot is the inferred maximum rotation speed 
and a z is the line-of-sight velocity dispersion averaged roughly over an Arecibo beam; see 
Paper I for discussion. In the case of Arecibo or other single beam maps, we determined 
these quantities as in Paper I; for synthesis array maps we determined the quantities from 
the rotation curves provided by the authors, or from galaxy-wide (or single-beam) H I 
profiles in a manner similar to that employed in Paper I. Since inclination corrections to 
the rotation speed become quite uncertain at small inclination angles, and self-extinction 
corrections to the blue luminosity become uncertain for nearly face-on galaxies, we will 
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mainly restrict the sample to 45° < i < 85° for types Sb through Sm. Im and BCD galaxies 
of all inclinations are retained in the sample since it is not clear how inclinations should be 
measured for these galaxies in any case. 

The diameter of the outermost recorded H I intensity contour is taken for D H ,max- For 
our Arecibo maps this point is assumed to be about half a beam-width in from the center of 
the outermost observed beam as explained in Paper I, and we defined D H ,max similarly for 
other single-beam maps. In the case of synthesis array maps for which a contour map was 
displayed, we simply measured the diameter of the lowest-level contour shown. We defined 
Dn, e , the diameter at 1/e of the central flux, by fitting a flat-topped exponential model 
to the observed beam fluxes as detailed in Paper I for our Arecibo maps. For synthesis 
array maps this quantity could be measured directly from the displayed contour maps by 
identifying the contour at a level 1/e times the central contour level. Single-beam maps 
by other authors gave us more difficulty in defining Dn, e - In some cases, e.g., van Zee, 
Haynes & Giovanelli (1995), a diameter defined in a different way turned out to be close to 
our Du e for several galaxies in common between us and those authors, and so we simply 
adopted their definition as equivalent to ours. For the several galaxies in the sample for 
which maps did not extend to the sensitivity limit we have determined Dn^max from the 
regression of Dn,max vs. Dn, e - In this paper we shall use Dn^max rather than Dn, e since 
the former is more directly related to the total mass of the galaxy. Unfortunately not 
all galaxies are measured to the same sensitivity limit; in particular, many of the early 
synthesis array maps do not reach below a sensitivity of 10 20 atoms cur 2 . This introduces 
considerable systematic error into the measurement of Dn,max in the sense that a number 
of the values are too low. 

3. Selection Biases, Extensive Variables and Comparison of Samples 

3.1. Selection bias and hydrogen radius 

The selection criteria stated above introduce a selection bias in some (but not all) 
properties. There is no bias against low optical surface brightness (if one discounts 
objects with such extremely small El that special detection techniques are required), and 
our sample ranges over a factor of ~ 10 3 in E^. There is considerable bias against unusually 
small hydrogen mass Mh for galaxies selected for hydrogen mapping. There is no direct 
bias against galaxies with very small hydrogen radius Rh, i.e., a galaxy is included (as 
long as mapping was attempted) even if it has only an upper limit for R H . However, since 
hydrogen surface density E# oc Mh/ Rh 2 tends to lie in a narrow range, the absence of very 
small M H indirectly discriminates against very small Rh- The M H bias has two effects: 
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(i) galaxies of any luminosity L, which are greatly gas-deficient due to ram-pressure or 
tidal stripping, are likely to be missing, and (ii) the lower end of the normal distribution of 
Mh/L is depleted for small L but not for large L (for indirect estimates of this effect see 
Sects. 4.2 and 5). 

As discussed in the preceding section, we have attempted to obtain a value for the 
"isophotal H I radius," Ru,max = \D H ^ max for each galaxy in our own sample from Paper 
I and in the larger sample from the literature in Sect. 2 of this paper. For a few galaxies 
in the combined sample we have only an upper limit to RH,max and in those cases we have 
taken the measurement to be half the formal upper limit. In some cases this will be an 
overestimate of the true radius. On the other hand, for a few galaxies in the sample drawn 
from the literature the adopted H I radius may be too small as discussed in Sect. 2. For 
all statistical purposes we shall omit the four "special" galaxies of Paper I along with HI 
1225+01 (Giovanelli & Haynes 1989) which has a uniquely large H I to optical radius ratio 
as we shall document below. On the other hand, our omission of galaxies with particularly 
large Rh is offset by the indirect bias against small Rh mentioned above, i.e., the tails on 
both sides of the distribution in R H /R opt are depressed somewhat. 

The optical isophotal radius R25 is known quite reliably for all the galaxies in our two 
samples, and Figure 1 displays Ru^ax vs. R25 for each galaxy in the full sample, (a) In the 
Figure, galaxies from Paper I are shown with larger symbols than those from the literature. 
The inclusion of spirals in the sample flattens the slope of the best-fit correlation slightly 
(although it is still within the uncertainty of the slope for the sample of Paper I alone), 
but there is no evident systematic difference between the radii presented in Paper I and 
those obtained from the literature, (b) As discussed in Paper I, many of the galaxies in our 
sample are drawn from the Virgo Cluster Catalog (Binggeli, Sandage & Tammann 1985) 
and some of those are within the core of the cluster, < 5° from its center. Similarly, some 
of the mapped spirals drawn from the literature are in the Virgo cluster core. While it is 
well-documented that some cluster spirals (Haynes, Giovanelli & Chincarini 1984) and 
irregulars (Hoffman, Helou & Salpeter 1988) are highly deficient in H I, such galaxies are 
not favored by our selection criteria and any effect of stripping is lost in the overall scatter 
of the points in Fig. 1. (c) Fitting the 116 spirals separately from the 109 dwarfs, we obtain 
(using bisector regressions as discussed below) for spirals Ru,max = (2.02 ± .24)i?25 980± '° 51 ' ) 
and for dwarfs R H ,max = (2.72 ± .ll)R^ 039±mi) . The powers are identical, within the 
uncertainties, but the fact that the scale factor is larger for dwarfs than for large spirals (by 
a factor of 1.4) is presumably due to two causes: (i) R H /R opt is overestimated, because of 
the selection bias, for small L (dwarfs) but not for large L (spirals), and (ii) galaxies with 
an irregular morphology may have less efficient star formation and may have retained more 
of their gas (see also Sect. 4.2). Linear regression on the logarithms of the radii (using the 
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bisector of the direct and inverse regressions; see Sect. 4 and Isobe et al. 1990) gives the 
correlation for spirals and dwarfs together: 

R / T3 \ -0.110±.031 

= (2.34 ± .14) 2 f (1) 
R 25 1 ; V 4 - 25k P c / 

where the normalization length, 4.25 kpc, is a typical value for i? 2 5 corresponding to 
L B = 1O 9 L . 



3.2. Definitions of extensive variables 

To convert angular radii and magnitudes to physical units, we have taken distances 
given by a Virgocentric infall model with asymptotic H = 74 km s~ 1 /Mpc and a Local 
Group deviation from Hubble flow of 273 km s _1 toward Virgo, as in Paper I. A number 
of the nearby dwarfs (especially those in the Local Group) have distance estimates from 
primary or secondary distance indicators; in those cases we have adopted the distances used 
in the sources of the H I mapping with a few exceptions noted in Table 3 of Paper I. The 
relative tightness of the correlation and its small deviation from linearity suggest that the 
hybrid surface densities (H I mass over optical radius squared, e.g.) used by various authors 
differ by only a scaling factor from quantities derived from H I measurements alone (with a 
few significant exceptions as indicated by the solid symbols in Fig. 1). In later sections we 
prefer to use a single radius for correlations against other variables. RH,max is in principle 
more appropriate since it gives the outermost point where we have velocity measurements, 
but i? 25 is measured more reliably. As a compromise, we shall use the geometric mean 
radius R gm between RH, ma x and 2.34i? 2 5, so that R gm is close to RH, ma x in most cases: 



Rgm - (2 .34^25)° l5 R°H,max- ( 2 ) 

We shall find it convenient to normalize the various extensive quantities, namely 
luminosity, radius, velocity and H I mass, by their values appropriate to a luminosity of 
10 9 Lq. Thus we define 



R 9m V c M H 



m H = r ^ Z m«u. • ( 3 ) 



1O 9 L ' ' 12.30kpc' 80.51 km s" 1 ' u 5.76 x 1O 8 M 

For the indicative dynamical mass Md yn = V c 2 R gm /G we then have 
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m dyn = v 2 r = M dyn /1.212xl0 10 M Q . (4) 

Although rridyn is not an independent observable, but simply constructed from v and r, we 
shall treat it separately so that we have five normalized extensive variables £, r, v, ran and 
mdyn- There are then 10 pairs of variables, and we display just 6 of the 10 in Figs. 2 and 3. 

We will discuss the correlation statistics in detail in Sect. 4, but first we note how the 
galaxies from Paper I compare to the larger sample. In each of the panels of Figs. 2 and 3, 
the Paper I data are shown with a larger symbol than data extracted from the literature. 
The powers in the correlations above are all slightly smaller than those given in Paper I for 
our data alone, but always within the uncertainties (3a rule) in the powers for the smaller 
dataset. In most cases the flattening of the regression lines seems to be due mainly to the 
inclusion of the larger spirals in the sample, but in the case of log v vs. log £ there is an 
apparent displacement of the points from Paper I with respect to those from the literature, 
in the sense that our sample has smaller rotation speed at a given luminosity. This is most 
likely due to the selection by some authors (especially van Zee, Haynes & Giovanelli 1995) 
of low surface brightness objects for mapping: In terms of surface brightness defined 
below (Eqn. 9) their 11 galaxies (which survive our inclination cuts) have a mean of logS^ 
equal to (—0.81 ± .19), compared with (—0.22 ± .05) for the 58 surviving galaxies from 
Paper I. 

4. Correlation Statistics for Five Extensive Variables 

4.1. The combined sample 

In the previous section we defined four physically independent variables, £, r, v, and 
m H , which are "extensive" in the sense that they are related to overall size and mass. 
One complication is that the measurement errors are not quite independent since £ and 
itlh are proportional to d 2 and r is proportional to d, where d is the adopted distance 
which is sometimes controversial. A second complication is related to m^ yn , defined in Eqn. 
(4), which is physically even more important as an extensive variable since total mass is 
less likely to fluctuate with time than £ does because of starbursts. The mass md yn can 
be treated as an independent variable for correlation purposes, but the results may be 
inconsistent with the fact that m^ n is also the product of v 2 and r. 

Consider first the choices we would have if md yn were truly independent. Physically, we 
have to consider all five extensive variables on an equal footing a priori, so in each of the 
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10 pairs of variables we should choose one of the forms of correlation analysis that treats 
the variables symmetrically. Variations in each variable are large, so that statistical results 
can depend strongly on whether the physical variables themselves or their logarithms are 
used. Especially because of measuring errors, the logarithms of the variables are closer to 
being normally distributed than the variables themselves. As in most previous analyses, we 
shall carry out all linear regressions using log 10 of the variables £, r, v, m H , and m dyn , so 
that raising a physical variable to a power p merely rescales the logarithm by a multiplying 
factor p (see Fig. 4 and Sect. 4b, however). Since we are interested in the fundamental 
relationship between the quantities rather than predictions for one quantity given a value 
for another, we seek a regression that treats errors in the two variables symmetrically. 
Isobe et al. (1990 — hereafter IFAB) have discussed in detail three different symmetric 
definitions of linear regression, have shown that they usually give different values for the 
slope f3 of a line y = a + f3x and have demonstrated that the "bisector method" (the chosen 
regression being the bisector of the direct and inverse ordinary least-squares regressions) 
is usually preferred. This (logarithmic) bisector regression for the other four extensive 
variables vs. t gives 

r = (1.002 ± .033)£ a382± - 013 , v = (1.001 ± .027)£ a276± - 014 , \ 

m H = (1.028 ± .075)£ a759± - 033 , m dyn = (0.999 ± .064)£°- 859± - 027 . J { } 

If the five variables were truly independent, a convenient way to proceed would be 
to keep the first power of one of the five physical variables and then to raise each of the 
other four to a particular power so that the rms deviation of the log 10 from its mean, a m , 
is the same for all five. We choose to keep the first power of £, for which a rn = 0.984 
in log 10 with the mean corresponding to I = 0.86 x 10 9 Lq (That a m is close to unity 
is purely a coincidence). The powers that give exactly a rn = 0.984 for the other four 
variables would be 2.639 for r, 3.683 for v, 1.322 for m H , and 1.165 for m dyn , and for those 
choices we would have unit slope for all 10 linear regressions no matter which of the three 
symmetric definitions we employ. This would be convenient, but leads to an appreciable 
inconsistency in m dyn = v 2 r, which is an example of the "non-commutativity" discussed 
in the Introduction: If r 2 639 and -y 3 - 683 were exactly equal to £, then v 2 = £°- 543 ; r = £°- 379 j 
and (i^r) 1 ' 165 = £ 1 - 074 . We thus have a 7.4% discrepancy in slope with the direct correlation 
between m dy ® 5 and £. 

Because of this discrepancy we cannot derive unique correlations but choose to adopt 
compromise powers which decrease the discrepancy slightly, although slopes are no longer 
unity. We choose 
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£, r 2 - 68 , v^, m]f 5 , m d f n (6) 

and use the bisector regression between all 10 pairs. The slopes all lie between 0.97 and 
1.03, and the formal statistical errors of the slopes (as defined in IFAB) lie between 3% and 
10%. Note that systematic errors in the slopes could well be much larger. The discrepancy 
in slope for v 2 x r vs. £ where v 2 and r are taken from the first three entries in Eqn. (5) 
against the direct m dyn vs. £ slope is 5.48%. The zero-points in the 10 regressions are 
within ±0.01 (in log 10 ) when the normalizations in Eqn. (3) are used, but systematic errors 
are likely to be much larger (although hopefully much smaller than a m 0.98). 

Although the slopes of the 10 regressions are all close to unity, the tightness of the 
correlations varies appreciably. Two measures of this tightness are displayed in Table 2: The 
numbers in brackets are the correlation coefficients (the geometric mean of the slopes of the 
direct and inverse ordinary least-squares regressions) while the numbers without brackets 
give the rms deviation a reg of either variable from the regression line. The various a reg are 
to be compared with a m ~ 0.98. The tight correlation between m^yn and v is mainly due to 
their interdependence {mdyn = v 2 r). On the other hand, the surprisingly tight correlation 
between uih and r 2 ' 68 is only partly due to the fact that errors in distance d almost cancel 
(rriH oc d 2 , compared with d 2S8 for r 2 68 ); the physical correlation must also be quite tight. 
The poorest correlation is seen between m]f 5 and t> 3 - 73 , but here we cannot tell which of two 
possible causes for the poor correlation is dominant: (1) Distance errors matter, since one 
variable is independent of d whereas mjf 5 oc d 2 ' 70 (and similarly, inclination errors affect v 
but not rnn), and (2) the physical correlation between mu and v might itself be poor. 

4.2. Possible nonlinearities or breaks between dwarfs and spirals 

The scatter in most variables is larger for the dwarf galaxies than for the regular 
spirals, partly due to genuinely larger fluctuations and partly due to warping and other 
asymmetry in the disks giving larger errors in inclination and hence in v. Linear regressions 
carried out separately for the dwarf irregulars (Hubble types Sdm and later) and for the 
regular late-type spirals (Sb to Sd) have larger statistical errors (especially for the dwarfs) 
than for the combined sample. We have nevertheless carried out these separate regressions 
for the 10 pairs of (log 10 of) the 5 variables in Eqn. (6), to look for any possible differences 
in slope and/or zero-point between dwarfs and spirals. Some of the results are displayed in 
the dashed lines in Figs. 1 and 2. Although there is some overlap in luminosity between 
types Sd and Sdm, results would have been very similar if we had separated £ < 1 from 
£ > 1 rather than by morphological type. 
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For the correlations involving H I, there seems to be a parallel displacement between 
the dwarfs and spirals, with the dwarfs being more H I-rich: For the hydrogen radius 
vs. optical radius relation in Fig. 1, in the vicinity of i? 2 5 ~ 3 kpc, the dashed lines 
give log 10 Rn,max larger for the dwarfs by (0.158 ± .074) than for the spirals. Similarly in 
Fig. 2a where 1.351og 10 m# is plotted vs. log 10 £, the dwarfs lie higher than the spirals 
by a zero-point difference of (0.71 ± .19), corresponding to a factor 3.4 in H I mass. The 
difference in log 10 RH,max corresponds to a factor 2.1 in (RH,max I R-25) 2 ■ The slopes of the 
pair of dashed lines are almost equal to each other (1.30 ± .09 for spirals, 1.22 ± .13 for 
dwarfs) but, because of the displacement, are both larger than for the combined single line 
(1 ± .03). The difference by a factor of 3.4 between dwarfs and spirals (or between faint and 
bright galaxies) is statistically significant only in a formal sense, and it is smaller than the 
variance in either group. Furthermore, this difference is reduced (although only slightly) if 
one corrects for the selection bias, discussed above, which has caused some galaxies below 
the regression line on the left hand side of Fig. 2a to be omitted. A downward correction 
factor of the convenient form (^/100) c (with positive c) should probably be applied to m# 
from our figures, but two arguments show that c is small: (i) Stavely-Smith et al. (1992) 
have carried out a dwarf-spiral comparison for m H vs. £ (without mapping) which had no 
selection bias on inn, and their power law for m# as a function of £ is very similar to that 
in our Fig. 2a with no c-correction; and (it) the selection bias affects only the low end of 
the distribution in mjf/f, but not the upper envelope of the scatter diagram in Fig. 2a. 
The slope of this upper envelope is not very different from our solid line, and we feel that c 
in the multiplicative correction factor (£/100) c for itlh could not be much larger than about 
+0.07. Our main physical conclusions are then that (a) there is no strong, sudden break 
between dwarfs and spirals, but (b) the ratio m H /£ decreases with increasing luminosity £, 
perhaps as steeply as £~ - 27 but probably more like £~ - 20 . 

Another correlation of interest is the Tully-Fisher (TF) relation, 3.73 log t> vs. log£, 
shown in Fig. 3a, where we again fit bisector regressions for spirals and dwarfs separately. 
A break in the zero-point would be expected if there were a physical distinction between 
dwarfs and spirals, such as massive neutrinos contributing to the dark halos of spirals but 
not to dwarfs. A mere change in slope would be expected if the TF relation is continuous 
but has curvature, as has been reported before. For spirals we obtain a slope of (1.31 ± .14), 
about 3(T steeper than that for dwarfs, (0.88 ± .10). The zero-points are almost the same, 
(—0.21 ± .14) for spirals and (—0.10 ± .11) for dwarfs, suggesting continuity plus curvature 
(a continuous change in slope). Although the errors in slope are large for our dwarfs and 
spirals separately, the significance of the change in slope is strengthened by comparing our 
combined sample (for which the error is much smaller, but dwarfs and irregular galaxies 
are still emphasized) with previous work: Writing £ oc v s , our total sample gives a slope 
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s = (3.73 ± .19), while Pierce & Tully (1988) find s = (2.74 ± .10) for regular (bright) 
spirals alone. For practical application of the TF relation, some authors use the direct 
regression for regular spirals alone, which gives an even smaller value of s ~ 2.5. However, 
for a physical understanding of the Hubble sequence, our use of symmetric regression and 
inclusion of dwarfs is more appropriate, and our averaged slope of s = 3.73 is closer to the 
original TF suggestion of s = 4. 

The regression of logf vs. logr in Fig. 3b also shows some curvature, so that the 
luminosity-radius relation in Fig. 2c has almost constant slope throughout, £ oc r 2M or r 2 7 
(the small displacement is probably not significant). For the three log-log relations shown 
in Fig. 3, the curvature can be represented by a two-component fitting formula, but the 
slopes and ratios of coefficients are very uncertain. The dashed curves in Fig. 3 show the 
following particularly simple choices: 

v 5 « 0A£ + 0.2£ 2 , (7) 
v 2 « 0.5r + 0.4r 2 - 5 , (8) 

m dyn w 0.4£ - 6 + OAf- 2 . (9) 

The use of logarithms in linear regressions is convenient and almost necessary if 
variables have a log-normal distribution. However, the implied "logarithmic averaging" may 
not be the physically meaningful procedure. In Newtonian mechanics the Virial Theorem 
requires that velocity-related averages over an orbit be carried out for v 2 rather than for 
any other power or for a logarithm. We illustrate this with the logarithmic v vs. £ relation 
in Fig. 3a for the dwarfs, where the dashed curve approximates t/ 4 - 22± - 48 ) = (0.76 ± .21)£ 
at the dwarf end. In Fig. 4 we have replotted the same data for dwarfs on linear scales for 
v 2 and £ 1/2 and show the bisector regression line v 2 = (-0.039 ± .048) + (1.02 ± A2)£ 1/2 . 
In principle, averaging v 2 over a log-normal distribution centered on the zero-point of 
the 3.73 log v - \og£ regression could give (v 2 ) significantly larger than the square of v 
corresponding to the zero-point itself, but in fact with the variance 0.652 in 3.73 logv we 
get (v 2 ) = 1.22, almost the same as the value 0.98 from the regression of v 2 vs. £ x l 2 . 



5. Three Surface Densities 



Three surface densities are of physical interest, namely Ex = Lb/ ^R 2 gm (related to 
the average surface brightness in the blue), T, H = M H / AnR 2 gm1 and E rf?/n = M dyn /4:TTR gm . 
These are related to the extensive variables in Eqn. (6) by 

Ex _ I S g _ m H ^dyn _ V 2 / lf y, 

1.23 L pc- 2 ~ r 2 ' 0.709 M pc- 2 ~ r 2 ' 14.9 M p C - 2 ~ r ' [ ' 

We should emphasize that the radius r, which we use in Eqn. (10) and in the rest of 
this paper, is mostly an "isophotal" or "isosurface-density" radius. If we could use some 
"scale-length" radius r s i instead (unfortunately unavailable where the mapping is coarsely 
resolved), the variation in Ex and would be greater, since low / high central surface 
brightness leads to a smaller / larger ratio of r/r s i (although only by a logarithmic factor 
for an exponential disk). Central values of surface brightness and surface density would 
presumably also have larger variations. 

If the correlation coefficients of the pairs of extensive variables were unity (i.e., straight 
lines on log-log plots and a reg = 0), and if the md yn inconsistency discussed in Sect. 4 
were absent, Eqn. (6) would give uniquely the dependence of the surface densities on 
the extensive variables. In reality, the correlations are not very tight, and there is some 
ambiguity and some curvature. To illustrate the uncertainties, consider first only the first 
four entries in Eqn. (6) and pretend that they represent exact equalities. This would give 



E L oc r - 68 oc £ a25 oc v°- m oc m° H 3 \ 

-0.015 „ fi-0.006 „ „, -0.021 ^ ^,-0.007 



Tih K r -0.015 K £-0.006 K ^-0.021 K m -U 

^ dyn oc r - 44 oc £ - 16 oc v 0M oc m° H 2 



;n) 



Instead, we can obtain the "direct logarithmic regression" for each of the three surface 
densities separately against each of the five extensive variables in Eqn. (6). The five 
alternative regression slopes (and their average) are shown in Table 3 for each of the three 
surface densities, together with the formal statistical error. The alternative slopes can differ 
by more than the formal errors, which illustrates a systematic artifact for poorly correlated 
quantities: A ratio like Ex, oc Ijr 2 has a larger slope in its regression against I and a smaller 
slope against r 2 ' 68 than in Eqn. (11), whereas the regression against a "neutral" extensive 
variable (such as v 3,73 or £r 2 ) is probably more reliable. 

To summarize our data so far on the variation of the surface densities along the mass 
sequence: The blue surface brightness definitely increases steadily (without curvature) along 
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the sequence, approximately as T, L oc (£r 2 ) e , with e ~ (0.13 to 0.16).. For our combined 
sampled, £# varies very little along the sequence, but we saw in Sect. 4.2 that a small 
correction in the form of a multiplying factor ~ (£/100) 07 should probably be applied 
to itlh and hence to E#) because of the bias against faint galaxies with small m#. That 
there is a correction to but that it is small, can also be seen indirectly from the scatter 
diagram in Fig. 5a of £# vs. £ l for our sample: The selection bias applies only to the lower 
portion of the diagram, not to the upper envelope which has only a small positive slope. 
Some of the positive correlation must be an artifact due to measuring errors in r, which 
enters both T, H and T, L to the same power. Thus, even the corrected T, H should increase 
little along the sequence. The variations in the scatter diagram of logS L vs. logE^ 
(Fig. 5c) combine systematic variations along the sequence with individual variations 
"perpendicular to the sequence." If one attempts to describe the "parallel" variation (along 
the sequence) as oc (Ed 2/n ) p n, then Eqn. (10) and Table 3 give p\\ ~ 1.6 or 1.8. However, 
the inclusion of curvature in Eqn. (8) gives 



H dyn oc (1 + 0.8r L5 ) (12) 

for the variation along the sequence, i.e., T: dyn hardly varies at the low end of the sequence 
(whereas S L does), and p\\ is not a very meaningful parameter. 

To investigate the variation of against T, dyn "perpendicular to the sequence," we 
eliminate the variations along the sequence by using slightly altered definitions: We choose 

yl _ El _ l°' 84 yi _ Z^dyn _ 

L ~ (£ r 2)0.16 - r 2.32' dyn- 1+ ~ r + Q g r 2.5 ^ > 

so that E' L and T! dyn have little systematic variation along the sequence. The scatter 
diagram is shown in Fig. 5d and the bisector regression line gives 



E' L oc(E^ n r, p± = 1.118 ±.085 (14) 

Note that a similar slope would be obtained if there were little true correlation between Y,' L 
and T>' dyn but large measuring errors in r. The true value of p± is therefore more uncertain 
than indicated by the formal error in Eqn. (14). 



6. Conclusion 
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6.1. Summary. 

We have analyzed the statistical relations of five extensive (size-related) observables 
for a large set of galaxies, with a broad range of absolute luminosity and surface brightness. 
Only galaxies with H I mapping were included, early-type spirals (SO, Sa, Sab) and 
ellipticals were excluded along with spirals thought to be members of tidally interacting 
pairs or groups. Faint dwarf irregular galaxies are more prominent in our sample than they 
would be in a magnitude-limited catalog. These five observables are t (a normalized blue 
luminosity Lb), r (a normalized "geometric mean" radius R gm , formed from an isophotal or 
isosurface-density H I radius close to that of the outermost detectable H I and the optical 
radius R25), v (a normalized velocity profile width V c which incorporates both rotation 
and random motion), uih (the normalized total mass Mjj of neutral atomic hydrogen), 
and m dyn (the normalized indicative total dynamic mass out to radius R gm , defined as 
M dyn = V 2 R grn /G). These quantities are normalized by 10 9 Lq, 12.30 kpc, 80.51 km s _1 , 
5.76 x 10 8 Mq, and 1.212 x 10 10 Mq, values typical of the region in which dwarf irregular 
and spiral galaxies overlap. 

For our complete sample we analyzed all Hubble types Sb through Im, which gives 
a particularly large dynamic range of the five extensive observables, which peak around 
Sbc (Roberts & Haynes 1994). Of the 225 galaxies in the complete sample, four were 
anomolously H I-rich, and these were omitted in all statistical analyses. Because of a 
selection bias for H I-mapping, extremely H I-poor galaxies are missing from the sample. 
As most previous authors have done, we correlated logarithms of the variables rather than 
the variables themselves. Since the five variables are on an equal footing physically a 
priori, we followed Isobe et al. (1990, IFAB) in using the "bisector" prescription for a 
symmetric form of a linear regression line for each of the 10 pairs of 5 variables. Because of 
the large scatter in the data, correlations suffer some ambiguity from "non-commutativity," 
especially for the product m^ n = v 2 r. Nevertheless, we find that the regression slopes 
for all 10 regressions between the log 10 of £, r 268 , t> 3 - 73 , mjf 5 , and m dy ® lie between 0.97 
and 1.03, with statistical errors <^ 10%. Because of the non-commutativity for ratios, we 
again have some ambiguities as to the variation along the size/mass/luminosity sequence 
for the surface densities = £/r 2 for blue luminosity, S# = m#/r 2 for neutral hydrogen 
mass, and T, dyn = v 2 /r for indicative dynamic mass. The extent of the ambiguity and one 
definition of an average slope for each pair are given in Table 3. We have also attempted to 
characterize two kinds of variations of with Y> dyn by first writing oc (S ( i J/n ) p n for the 
average changes along the mass sequence. These two surface densities were then corrected 
for the average systematic variation along the sequence to obtain new variables T! L and 
T/ d . Then we fitted a regression line of the form T! L oc (^' dyn ) p± to the corrected variables 
to represent the variation "perpendicular to the mean sequence." The regressions gave 
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p\\ ~ 1.6 or 1.8 and p± ~ 1.15, but these present values are unrealistic for reasons given 
below. 

We have looked for curvature in the logarithmic regressions of the various extensive 
variables against each other, and we also investigated whether or not there are any 
"breaks" in the correlations for "dwarfs," defined as Sdm through Im (including BCD), 
irrespective of size or luminosity, vs. those for "spirals," types Sb through Sd. For the 
variables not involving H I, we find no significant breaks between Hubble types, but the 
large range of extensive variables enabled us to find curvature in some of the regressions: 
For the Tully-Fisher relation £ oc v s , the bright spirals (which are preferred for practical 
application) give s ~ 2.5 or 2.7, our total sample gives s ~ 3.7, and the fainter galaxies 
give s slightly larger than 4, the power originally suggested by Tully and Fisher (1977). 
There is little curvature in the relation £ oc r 2 - 7 , so that the surface brightness in the blue, 
El = £/r 2 , increases steadily along the mass sequence (El in the red or infrared increases 
even more rapidly). There is curvature in the logw-logr relation, so that the surface 
density of dynamic mass is almost constant at the low end of the mass/radius sequence 
and increases appreciably along the upper end of the sequence. This variation, roughly 
^dyn — v 2 /r oc (1 + r ' 8 ), coupled with the steady progression of El, shows that a relation 
of the form El oc (S^„) p n is not very meaningful. The concept of the "perpendicular 
parameter" p± in E' L oc (E' dyn ) p± is meaningful, but out present value for it (~ 1.15) may 
be quite uncertain because of a coincidence involving powers of the galactic radius r, which 
may have large measuring errors at the moment: The definitions of E' L and E' dyn are such 
that particularly large errors, in the absence of strong real correlation, could mimic our 
results. 

Only for correlations involving H I content are there tentative indications of "breaks" 
between dwarfs and spirals in the sense that the Rn-R pt relation lies higher by a factor 
of ~ 1.3, and the m H -£ relation by a factor of ~ 3, for dwarfs than for spirals. These 
"displacements" near I ~ 1 are smaller than the variances and it is not clear if the 
effect is real. The surface density of H I mass, in our sample varies little along the 
luminosity/mass/radius sequence, and the individual variations from En ~ 0.7 Mq pc~ 2 
are fairly small. Particularly small values of mn/£ and of Eh tend to be absent (especially 
for dwarfs) because of a selection bias, and we have estimated that mn/£ and Eh should 
be corrected downwards from our results by a multiplying factor of order (Lb/IC^Lq) - 07 . 
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6.2. Comparison with previous analyses. 

Four recent papers have carried out related correlation studies with slightly different 
emphases: (i) Roberts & Haynes (1994) mainly investigated variations with Hubble type 
rather than with luminosity, (ii) Gavazzi (1993) concentrated on variation with luminosity 
as we do, but he included SO, Sa and Sab galaxies and did not have much data on very 
faint dwarfs or low surface brightness (LSB) galaxies. On the other hand, (Hi) Sprayberry 
et al. (1995) and (iv) Zwaan et al. (1995) present extensive data on LSB galaxies and 
compare with galaxies (of similar dynamic mass M dyn ) with normal surface brightness S^. 
While there is agreement on most trends, there is considerable controversy regarding the 
trends of (in the blue) and T, dyn along the size / mass / luminosity / Hubble sequences 
and also on their correlations relative to each other "perpendicular to the sequence" (i.e., 
variations for fixed m dyn or £). We show first that the apparent discrepancies are partly due 
to uncertainties brought about by some near cancellations for variations along the sequence: 
Assume that £ = v a and r = v 13 along the whole mean mass sequence. With Ti dyn = v 2 /r 
and S L = l/r 2 , we have 

X dyn = S L = n = ^l. (15) 

Our derived approximate values are a ~ 3.8 and [3 ~ 1.4, so that a — 2(5 is appreciably 
smaller than a and 2-/3 smaller than 2. This leads to a large uncertainty in p\\ even for 
constant a and (5. Furthermore, we saw that a and (3 vary along the sequence (with 2-/3 
essentially vanishing on the low end), so that p\\ is not well defined. 

Consider next the "variations perpendicular to the sequence," i.e., changes in E^ n at 
fixed m dyn = v 2 r, so that v 4 oc r~ 2 oc T,^ yn . Assume that S L oc (S d?/ri ) p± for variation at 
constant m dyn with some value for the constant p±, so that 

£ OC (Edyn) P± ~ 1 , i/v a = (^ d yn) {P± - 1] - a/4 . (16) 

For the special case of p± — 1, the luminosity / mass ratio would then not depend on T> dyn 
variations. For the special case of p± = 1 + a/4 on the other hand, variations in Y, dyn away 
from the mean value along the sequence would not affect the Tully-Fisher v-£ relation at all. 

Both Zwaan et al. (1995) and Sprayberry et al. (1995) compared the Tully-Fisher 
(TF) relation between v and I for normal galaxies with that for LSB galaxies, which gives 
some information on Eqn. (15). However, there is some spread in the observations and the 
two papers emphasize different aspects: The former notes that the difference in the two 
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TF relations is small, so that p±_ should be close to (1 + a/4), which is about 1.9 (in the 
blue and using our bisector slopes). On the other hand, the latter paper notes that there 
is some difference in the TF relations — in the sense that "Malin I-like" objects with very 
large r and small v, small E^, have anomalously small £/v a . These few galaxies at least 
indicate that p± < 1.9 (If p± were as small as 1, then £/md yn would depend little on E^n, 
but £/v a would have a strong dependence.). Qualitatively, our own result of p± ~ 1.12 is 
intermediate between 1 and 1.9, as it should be, but the combined data of Zwaan et al. 
(1995) and Sprayberry et al. (1995) suggest a slightly larger value for p± than ours. We 
saw that the large present-day errors in radius make the value of p±_ uncertain. 

There is only an illusion of a discrepancy regarding the ratio £/m dyn , which Gavazzi 
(1993) has decreasing along the mass sequence whereas our averaged linear regression has 
m dyn- Our Eqn. (9) indicates "regression curvature" in the sense that £/m dyn 
increases along the sequence for small £ and then decreases slightly for large £ [Fig. 2d in 
Roberts & Haynes (1994) also indicates a peak in £/md yn at an intermediate Hubble type]. 

6.3. Discussion. 

An important negative result of our compilation is the absence of any sudden break 
between dwarf irregulars and regular spirals for variables not involving H I. In an overall 
dynamic sense, small dwarfs are thus a continuation of the spiral Hubble sequence. Only 
the fact that the Rn-Ropt and mn-£ relations seem to lie higher for dwarf irregulars than 
for spirals might indicate that gas depletion by star formation is reduced by morphological 
irregularities. However, even this fact may only be a manifestation of the tendency 
(discussed below) that hydrogen content decreases more slowly with decreasing mass than 
luminosity does. 

Of the three surface densities discussed above, only the blue surface brightness 
increases steadily and appreciably along the mass/radius/luminosity/Hubble sequence. A 
fourth density E* for total stellar mass is not measured explicitly for most of our sample, 
but should be correlated best with Ex, in the visible to near infrared. We do not have E^ 
at those wavelengths either, but L, and hence E L , is known to increase more rapidly along 
the sequence at longer wavelength than in the blue. Our uncorrected data has neutral 
hydrogen mass surface density E# almost constant along the sequence; the correction 
for selection bias would presumably have E# increase along the sequence, but only very 
slowly, and Ejx/E* definitely decreases along the sequence from dwarfs to giants. As a 
consequence, the total disk mass density (E# + E*) increases slowly (and is gas dominated) 
along the lowest portion of the sequence and increases rapidly along the middle and upper 
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portions of the sequence. The density E d?/ „ of total dynamic mass (mostly dark matter) 
is almost constant along the lower sequence and increases along the upper portion, so 
that the ratio (E# + E^/E^n probably increases quite slowly along the sequence. The 
small variation of E#, compared with the large variation of E*, fits in with the suggestion 
by Kennicutt (1989) that star formation proceeds rapidly when the gas surface density 
exceeds a threshold value: Rapid star formation on the upper sequence (where the initial 
gas density E# + E* is large) depleted the gas rapidly, but the rapid depletion did not occur 
on the lower sequence so that the present-day E# is almost constant along the sequence. 

The small variation of E# in our Figs. 5a and b stems in part from our radius definition 
and from a selection effect: As mentioned in Sect. 4, the use of isophotal radii r (instead 
of scale-length radii r s i) decreases variations in both E^ and E#, compared with central 
values (de Blok, van der Hulst & Zwaan 1995, private communication). Galaxies with very 
small E# are likely to have had low total m H and therefore not to have been chosen as 
candidates for mapping in H I, so the lowest portion of Figs 5a and 5b may be absent from 
our sample. However, this is not the case for upward fluctuations, and it is significant that 
only a few out of the 225 galaxies have an unusually large H I content for their position 
on the mass/luminosity sequence. This also meshes with the notion that an initially high 
gas surface density usually produces rapid star formation and eventually depresses E# to 
"typical values." 

We also have some qualitative information on the variation along the mass sequence 
of two different volume densities: (i) The volume density uh of H (and other gas) in the 
midplane of the galactic disk depends on E# x Ti dyn and inversely on the velocity dispersion 
A Vj_ normal to the plane. Although V c is smaller for dwarfs than for large spirals, AV± varies 
little, as does E#. Since T, dyn decreases slowly with decreasing size/mass/luminosity, n# 
also decreases slowly, (ii) For studying tidal effects on a galaxy, another quantity with the 
dimensions of a volume density is important, namely n tid = M dyn oc E^/ 'R gm . 

Using Eqn. (6), we see that this quantity is actually slightly larger for dwarfs than for 
spirals, varying roughly as n tid oc £~ a26 . Tidal effects should therefore be less pronounced 
on dwarfs than on spirals. (Hi) Although dwarfs are more "robust" from the point of view 
of (ii), galactic winds, "blow-outs," etc., depend on escape velocity [oc (T'dynRgm) 1 ^ 2 ], and 
dwarfs are less robust from this point of view. 

Since the extensive variables peak around Sbc, it would be of interest to explore 
the differences in correlations of the extensive variables and surface densities for three 
morphological groupings: Sdm through Im (including BCD), Sbc through Sd, and SO 
through Sb. There are at present too few detailed H I mappings of early-type spirals to 
allow this comparison between the three groupings, but mapping is proceeding at a great 
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rate. It would also be useful to make more detailed comparisons between galaxies with 
and without H I mapping to obtain p±_ more reliably. For discussions of the empirical i-v 
Tully-Fisher relation one wants to know if p± is always close to (1 + a/4); if so, that would 
be a partial explanation for the narrow width of the relation. 

We acknowledge fruitful discussions and correspondence with G. Bothun, J. Charlton, 
E. Corbelli, E. de Blok, J. Dickey, G. Gavazzi, G. Helou, S. McGaugh, M. Roberts, E. 
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Fig. 1. — Plot of radius RH,max to the outermost measured neutral hydrogen vs. optical 
radius i? 25 at the 25 mag arcsec~ 2 isophote. Irregular galaxies from Paper I are shown as large 
exes with BCDs as large open squares and upper limits shown at one-half the formal limit 
as open triangles. Mapped irregulars from the literature are shown as small exes or squares, 
similarly. Members of interacting binary systems are shown as asterisks and excluded from 
all correlations. Mapped spiral galaxies from the literature are shown as small dots. Open 
circles mark irregular galaxies with high luminosity (> 10 10 Lq). Five special cases are 
indicated with filled symbols and are omitted from all correlations: DDO 154 by a circle, 
DDO 137 by a triangle, the NGC 4532 / DDO 137 complex by a hexagon, VCC 2062 by 
a diamond, and HI 1225+01 by a small square. The solid line is the bisector of the two 
ordinary least-squares regression lines (y vs. x and x vs. y). The two dashed lines are 
bisector regression lines for the spirals and dwarfs separately. 

Fig. 2. — Logarithmic scatter plots and regressions of pairs of extensive variables, each 
normalized to the value corresponding to a blue luminosity of 10 9 Lq as described in the text. 
The powers (multiplicative factors in the logarithms) are explained in Sect. 4a. Symbols 
are chosen as for Fig. 1. The solid line in each case is the bisector regression line, while 
dashed lines are for dwarfs alone or spirals alone, (a) Neutral hydrogen mass m# vs. blue 
luminosity £. (b) m# vs. indicative dynamical mass md yn - (c) Geometric mean radius r vs. 
t. 

Fig. 3. — Logarithmic scatter plots and regressions of pairs of extensive variables, each 
normalized to the value corresponding to a blue luminosity of 10 9 Lq as described in the text. 
The powers (multiplicative factors in the logarithms) are explained in Sect. 4a. Symbols 
are chosen as for Fig. 1. The dashed curve in each panel is an a smooth curve which 
approximately matches the regressions for dwarfs alone at one end and for spirals alone at 
the other, (a) Velocity profile half-width v vs. blue luminosity £. (b) v vs. geometric mean 
radius r. (c) Indicative dynamical mass m^ n = v 2 r vs. £. 

Fig. 4. — Scatter plot and regression of normalized velocity profile half-width squared vs. 
square-root of normalized blue luminosity for dwarfs alone on linear, rather than logarithmic, 
axes. The solid line is the bisector regression for the displayed points. 

Fig. 5.— Logarithmic scatter plots of neutral hydrogen surface density, E# = m#/r 2 , 
vs. (a) optical surface brightness, E^, = i/r 2 and (b) dynamical mass surface density, 
^dyn — v 2 /r, and of (c) optical surface brightness, Ex, = i/r 2 vs. dynamical mass surface 
density, E^ n = v 2 /r and (d) reduced optical surface brightness, E' L , vs. reduced dynamical 
mass surface density, E^ , as discussed in the text (Sect. 5). Symbols are defined as in Fig. 
1. 



